Hydrogenated grain boundaries in graphene 
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We have investigated by means of first principles calculations the structural and electronic prop- 
erties of hydrogenated graphene structures with distinct grain boundary defects. Our total energy 
results reveal that the adsorption of a single H is more stable at grain boundary defect. The elec- 
tronic structure of the grains boundaries upon hydrogen adsorption have been examined. Further 
total energy calculations indicate that the adsorption of two H on two neighbor carbons, forming a 
basic unit of graphane, is more stable at the defect region. Therefore, we expect that these extended 
defects would work as a nucleation region for the formation of a narrow graphane strip embedded 
in graphene region. 

PACS numbers: 73.22.-f, 73.22.Pr, 71.55.-i 



Physical and chemical properties of hydrogenated two- 
dimensional graphene-based (nano) structures have been 
the aim of recent researches. The adsorption of H atoms 
can modify the structural, electronic and magnetic prop- 
erties of the ideal graphene, which is a null-gap semi- 
conductor [l|. It is known that the adsorption of a sin- 
gle H at graphene opens an energy gap and introduces 
spin-polarized gap states l2|- On the other hand, a fully 
hydrogenated graphene [3] is an insulator material [4], 
known as graphane. Partially hydrogenated graphene 
structures, such as graphene/graphane superlattices [5], 
graphene nanoroads [6] and graphane nanoribbons [7] 
have been theoretically proposed, and can be produced 
by removing H atoms from graphane [H or adsorbing 
them on graphene [ol, [lo| . These structures can present a 
suitable band gap energy [ol, [lo| for pratical applicantions 
in the new nanoelectronic devices. 

Grain boundary (GB) defects have been observed in 
graphene in very recent experiments JT, Hj. Theoreti- 
cal works [l3| have also proposed structural models for 
GB defects, which resemble to those of Refs. [ll|, Il4|. 
Tight-binding calculations suggest that the H adsorp- 
tion is more stable at the GB defects than in perfect 
graphene [15]. However, there is no ab-initio investiga- 
tions on the properties of H adsorption at those observed 
GB defects, H/GB. Thus, in this work we have investi- 
gated, by means of first-principles calculations, the effect 
of H adsorption on the properties of graphene with grain 
boundary defects. Our results indicate that the adsorp- 
tion of a single H is more stable at GB defect. Also, the 
adsorption of two H on two neighbor carbons, forming 
a basic unit of graphane, is more stable at the defect 
region. Thus, we can infer that these extended defects 
would work as a nucleation region for the formation of 
a narrow graphane strip embedded in graphene region. 
The modifications in the electronic structure due to H 
adsorption at GB defects are also investigated. 

Spin-polarized density functional theory (DFT) [l6| 
calculations, within the generalized-gradient approxima- 
tion (GGA) [17], were performed using the SIESTA 



code [18], where core states were replaced by norm- 
conserving pseudopotentials [l9[ in the factorized 
form [20] . We have employed a double-zeta plus polariza- 
tion (DZP) basis set with an energy shift of 100 meV [2l|, 
and an energy cutoff of 200 Ry for the real-space mesh. 
10 (210) special k-points for geometry (band structure) 
calculations were used. All atomic positions were fully re- 
laxed until the residual forces were converged to within 
10 meV/A. We have employed a distance of 15 A be- 
tween graphene layers and supercells with (60) 60, 96, 
and 50 carbons atoms to represent the (perfect) defec- 
tive graphene structures. 

We first describe the atomic structure of GB defects 
in graphene. Here, we have considered three distinct 
structural models for the GB defects, shown in Fig. [TJ 
These extended defects can be built from a suitable reg- 
ular arrangement of 5-, 7- or 8-membered carbon rings 
along a line in graphene. The GB structure shown in 
Fig.[T](a), named as GB 1(5-7), is formed by an alignment 
of pentagon- heptagon pair separated by hexagons (l4| . 
Figure. [T] (b) shows the GB structure named as GB2(5- 
7), which is a regular sequence of pentagon- heptagon 
pair [l3|. In these structures, due to the periodicity of the 
supercell, we have considered two GBs in the unit cell. 
The third GB structure, shown in the Fig.[T](c), is named 
as GB(5-8) and consists of a sequence of pentagon-pairs 
and octagon along the graphene zig-zag direction [12j . 

Let us now investigate the energetic stability of H 
atoms adsorbed at graphene with GB defects, H/GB. 
We first consider a single H atom adsorbed on the top of 
several non-equivalent carbon atoms at the defect region. 
These adsorption sites are indicated by the letters A, B, 
C, D, E and F, as shown in Fig. [TJ Also, we have consid- 
ered the H adsorption at the graphene-like region (site 
G). To investigate that, we calculated the adsorption en- 
ergy, Eads^ which is defined by the following equation: 



Eads = E[R/GB] - E[GB] - E[R] + (5bsse (1) 

where E[}i/GB] (^[GB]) is the total energy of the fully 
relaxed hydrogenated (pristine) graphene structure with 
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FIG. 1: Optimized geometries of distinct grain boundary de- 
fects, labeled as: (a) GBl(5-7), (b) GB2(5-7), and (c) GB(5- 
8). The dashed- line box indicate their periodic unit cells. The 
unit cell parameters a and b, in A, are 24.5 and 6.6; 28.9 and 
9.0; and 27.3 and 5.0 for GBl(5-7), GB2(5-7), and GB(5-8), 
respectively. The adsorption sites are in gray. 



GB defects, and ^[H] is the total energy of the atomic 
hydrogen. The last term (^bsse) is used to correct the 
basis set superposition error [23]. A negative value for 
Eads indicates that the adsorption of H is an energeti- 
cally favourable process. To test our methodology, we 
have calculated the adsorption energy for a single H at 
perfect graphene. The calculated adsorption energy of 
-0.81 eV/H atom is in excelent agreement with other ab- 



initio calculations 



TABLE I: Calculated adsorption energies, in eV/H atom, for 
a single H atom adsorbed at graphene with grain boundary 
defects. The letters indicate the adsorption sites (see Fig. [T]). 

H/GBl(5-7) H/GB2(5-7) H/GB(5-8) 



G 


-0.83 


-0.93 


-0.78 


A 


-1.55 


-1.75 


-2.29 


B 


-1.39 


-1.54 


-1.89 


C 


-1.30 


-1.45 


-1.59 


D 


-1.26 


-1.24 


-1.43 


E 


-1.18 


-0.99 


-1.31 


F 


-0.96 







Our calculated adsorption energies, shown in the Ta- 
ble in indicate that Eads at the graphene-like region (site 
G) are within 0.1 eV to that at the perfect graphene 
sheet. Thus, the error due to the size cell effect is small 
compared to our calculated energy differences. It is inter- 
esting to note that the adsorption energy at the GB de- 
fect is negative and smaller than that at perfect graphene 
for all investigated adsorption sites. This indicates that 
the adsorption of a single H is energetically more sta- 



ble and likely to occur at the defective region. Indeed, 
we find that H/GB(5-8) is a quite stable structure, since 
the hydrogen adsorption energy (2.29 eV/H) is compa- 
rable with the binding energy of an isolated H2 molecule 
(2.16 eV/H). The flatness of pristine graphene struc- 
tures is lost due to the H adsorption. For H/GB2(5- 
7), the hydrogenated carbon and its nearest neighbors 
move upwards by 0.5 and 0.2 A, respectively, while the 
other atoms exhibit negligible displacements. Meanwhile, 
for the H/GBl(5-7) and H/GB(5-8) systems, the hydro- 
genated carbon moves upwards by 0.4 A relative to 
its nearest neighbors, and a significant displacement is 
viewed even for the fourth nearest atom of the hydro- 
genated carbon. At the equilibrium geometry, the C-H 
bond lenght is 1.13 A and the bond length of hydro- 
genated carbon and its nearest neighbors ranges from 
1.49 to 1.55 A, indicating a sp^ hybridization for the car- 
bon atom bonded to H. 

Focusing on the electronic properties of the H/GB 
systems, we performed band structure calculations and 
scannning tunneling microscopy (STM) simulations for 
the most stable configurations. Figures [2] (a)- (c) present 
the electronic band structure of the pristine GBl(5-7), 
GB2(5-7), and GB(5-8) systems, respectively. The va- 
lence and conduction bands exhibit a linear energy- 
momentum relation nearly of the Dirac point for wave 
vectors parallel to the GBs. Our STM simulations, for 
the occupied states within an energy window of 0.2 eV 
below the calculated Fermi level [Figs. [3] (a)-(c)], indi- 
cate the formation of bright lines parallel to the GBs, 
which is in accordance with the more reactive charac- 
ter of the GBs in comparison with the perfect graphene 
sheet. Those findings are in good agreement with pre- 
vious works 12-1^. Figures [2] (d)-(f) show the elec- 
tronic band structure of H/GBl(5-7), H/GB2(5-7), and 
H/GB(5-8) structures, respectively. Our results show 
that the single H adsorption has a strong effect on the 
electronic properties of GB defects. The systems remain 
metallic, however, the Dirac point disappears, and it is 
noticeable the formation of spin-polarized states crossing 
the Fermi level for the hydrogen adsorbed GB 1(5-7) and 
GB(5-8) systems. In contrast, although the structural 
similarity with the H/GBl(5-7) system, H/GB2(5-7) ex- 
hibits a spin paired state crossing the Fermi level. 

The STM images of the GBs are modified upon hy- 
drogen adsorption. Figures [3] (d)-(f) show alternating 
bright and dark spots along a direction parallel to the 
GBs. The presence of H adatoms somewhat localizes the 
electronic states, which is in accordance with the small 
energy dispersion of the electronic states (near the Fermi 
level) as depicted in Figs. [2] (d)-(f). In H/GBl(5-7) we 
find the formation of bright spots on the hydrogen ad- 
sorbed site as well as on the nearby C atoms. Fig. [3] (d). 
Meanwhile, the STM image of the H adatom becomes 
faded for H/GB2(5-7) and H/GB(5-8). In H/GB(5-8), 
Fig. [3] (f), we find a very weak spot lying on the H 
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FIG. 2: Calculated band structure for pristine and hydro- 
genated graphene with GB defects, a) and d) GB 1(5-7), b) 
and e), (c) and f) GB(5-8). Solid (dashed) lines indicate the 
spin-up (spin-down) components. 



adatom, while the nearest neighbor C atom (parallel to 
the GB) becomes brighter. Indeed, our calculated pro- 
jected density of states for H/GB(5-8) indicate that the 
electronic states along the C-H bond are mostly local- 
ized between 2 and 3 eV below the Fermi level, while in 
H/GBl(5-7), the electronic states of H adatoms (mostly) 
lie at Ef ± 0.3 eV. In this case, in contrast with the 
H/GB 1(5-7) system, the contribution of the H adatom 
in H/GB (5-8) is not relevant for the spin polarized states 
near the Fermi level. Fig. [2](f). The STM results for the 
H/GB2(5-7) system indicate that the electronic states are 
less localized. Figs. [3] (e). The hydrogen adsorbed sites 
as well as the nearest neighbor C atoms present almost 
the same brightness along the GB. 

Finally, we consider the adsorption of two H atoms on 
graphene structures. We focus on the most stable config- 
urations for the single H adsorption. Thus, we place the 
H atoms on sites A and A' for GBl(5-7) and GB(5-8), 
and sites A and B for GB2(5-7). The configuration with 
the H atoms sitting on the opposite sides of the sheet 
is energetically more favorable than that on the same 
side. Also, our calculated adsorption energy on the per- 
fect graphene, -1.7 eV/H atom, is in excellent agreement 
with other ab-initio calculations [24]. For the defective 
graphene region, we find adsorption energies of -1.86, - 
1.98 and -2.12 eV/H atom for GBl(5-7), GB2(5-7) and 
GB(5-8), respectively. Therefore, the adsorption of two H 
atoms is more stable at the defect region. The structural 




FIG. 3: (Color online) Simulated STM images for the pristine 
(left panels) and for hydrogenated (right panels) GB defects 
(whitin an energy interval of (Ef — 0.2) eV and Ef). The 
densities are in electrons /bohr^ . 

configuration of the hydrogenated carbons resembles the 
basic unit of graphane. The above results allow us to in- 
fer that the formation of graphane stripes, embedded in 
graphene, can be ruled by the presence of grain boundary 
defects. 

In conclusion, we have investigated the effect of adsorp- 
tion of H atoms in the electronic and structural proper- 
ties of graphene with (extended) grain boundary defects. 
Our ab-initio calculations indicate that the GB defect 
could work as nucleation region for graphane strip for- 
mation. The electronic structure of such hydrogenated 
defect is strongly affect by the adsorption of a single H 
atom. Spin polarized states are observed in H/GBl(5-7) 
and H/GB(5-8) 
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